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ABSTRACT 

Context. Infrared Dark Clouds (IRDCs), condensed regions of the ISM with high column densities, low temperatures and high masses, 

are suspected sites of star formation. Thousands of IRDCs have already been identified. To date, it has not been resolved whether 

IRDCs always show star formation activity and, if so, if massive star formation (> 8Mq) is the rule or the exception in IRDCs. 

Aims. Previous analysis of sub-millimeter cores in the cloud MSXDC G048. 65-00. 29 (G48.65) indicates embedded star formation 

activity. To characterize this activity in detail, mid-infrared photometry (3-30 yum) has been obtained with the Spitzer Space Telescope. 

This paper analyzes the point sources seen in the 24 jira band, combined with counterparts or upper limits at shorter and longer 

wavelengths. 

Methods. Data points in wavelength bands ranging from 1 /jm up to 850 pm (Spitzer IRAC, MIPS 24 and 70 pm, archival 2MASS 

data and sub-millimeter counterparts) are used to compare each 24 pm source to a set of Spectral Energy Distributions of Young 

Stellar Object (YSO) models. By assessing the models that fit the data, an attempt is made to identify YSOs as such and determine 

their evolutionary stages and stellar masses. 

Results. A total of 17 sources are investigated, 13 of which are classified as YSOs, primarily - but not exclusively - in an early 

embedded phase of star formation. The modeled masses of the central stellar objects range from sub-solar to ~8 Mq. Every YSO is at 

less than 1 pc projected distance from its nearest YSO neighbor. 

Conclusions. IRDC G48.65 is a region of active star formation. We find YSOs in various evolutionary phases, indicating that the star 

formation in this cloud is not an instantaneous process. The inferred masses of the central objects suggest that this IRDC hosts only 

low to intermediate mass YSOs and none with masses exceeding ~8 Mg. 

Key words, stars: pre-main sequence - stars: formation ~ ISM: clouds - ISM: individual objects: MSXDC G048. 65-00. 29 
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1. Introduction 

Infrared dark clouds (IRDCs), discovered independently by 
Perault et al. ( 1996[ l and Egan et al. ( 1998[ l as dark patches 
against the Galactic mid-infrared background, are suspected to 
be sites of clustered star fc)rmation (Rathborne et al.|2006;,Simon 
et al. 2006b ). Moreover, the suggestion is raised (see e.g. Menten 
etaI.|2005t|Beuther & Steinacker|2007 1) that IRDCs harbor mas- 
sive young stars. Morphologies of IRDCs range from globu- 
lar to more filamentary structures. One of the most important 



properties of IRDCs is their high column density (> 10 cm 



e.g. Carey et al.|2000 1. Dust absorbs mid-infrared radiation and 
makes the clouds stand out as dark features against the mid- 
infrared background. Temperatures in IRDCs are < 25 K and 
volurne deii sities can exceed 10^ cm"^ ( Egan et al.|l998 [Carey 
et al.||2000|l. Typical size and mass scales are ~5 pc and 10^- 



the absence of a bright mid-infrared background), using 2MASS 
(Two Micron All Sky Survey) color distributions of background 
stars, followed up by observations of molecular lines and Spitzer 
Space Telescope photometry. 

It has not yet been established whether all IRDCs show ac- 
tive star formation or that some may in fact harbor only starless 
cores. Moreover, the association of IRDCs to massive (> 8Mq) 
star formation in particular is still a matter of dispute. Studies of 
star forming regions in general (e.g. Indebetouw et al. 2007| l and 
IRDCs in particular (e.g.'Beuther & Steinacker|2007 1 prove that 
star formation activity in a wide stellar mass range and at various 
phases of evolution can be probed by Spitzer photometry. 

The IRDC under investigation in this paper is the cloud at 
Galactic coordinates (C, b) - (48°66, -0?30). It has been pre- 
viously observed by the Mid-course Space Expeiiment (MSX), 
the SC UBA instrum ent on the James Clerk Maxwell Telescope 



KTMo (Simon et al. 2006b I. These low initial temperatures, (JCMT, Ormeletal. 2005), and by the JCMT in CO, '^CO 



high molecular densities and total masses typical for IRDCs are 
precisely what make them suitable candidates for star forming 
regions. IRDCs are generally not quiescent: they harbor compact 
cores of sub-millimeter emission (e.g. Carey et al.|2000 Ormel 
et al.|2005[ [Rathborne eTaL][2005| ). IRDCs are found primarily 
in the inner Galaxy and close to the Galactic plane. Thousands 
are known in the first and fourth quadrants of the Galactic plane 
(Simon et al. 2006a i, i.e., the inner Galaxy. Frieswijk et al. 
( 2007| l recently identified the first IRDC in the outer Galaxy (in 



and HCO"^ (Shipman et al.|J20031|). In addition, it is covered by 



the Galactic Legacy Infrared Mid-Plane Survey Extraordinaire 
(GLIMPSE, [Benjamin et al.|2003| l. In the dark cloud catalog of 



Simon et al.| ( |2006a[ ) this cloud has the designation "MSXDC 
G048.65-00.29", hereafter simply "G48.65". Its distance is de- 
termined kinematically from molecular line data and an assumed 
Galactic rotation curve. G48.65 is found to be at a distance of 
-2.5 kpc ( [Ormel et aL[|2005[ [Simon et aT][2006b| l. Its distance 
to the Galactic Center is ~1 kpc and it is less than 20 pc away 
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from the mid-plane of the Galaxy. Its total mass is estimated at 
almost 600 M p, within a 2 pc area a nd the molecular (H ?) density 
is -10^ cm"^ ( jSimon et al.|2006b| . |Ormel et al.| ( [2005] l identified 
three distinct emission cores at 450 and 850 //m; their modeling 
indicates the presence of central luminosity sources on the order 
of lO^-lO"' Lq in at least two of these cores. 

The earliest phases of star formation show emission at 
>30 //m as they are embedded in an accreting envelope; the 
peak wavelength shifts through the mid- and near-infrared to 
the optical regime as stars evolve toward the main sequence. 
As mentioned above, sub-millimeter observations indicate on- 
going embedded star formation in two or three cores in G48.65 
(Ormel et al. 20050. This raises the question whether the star 
forming activity in G48.65 is limited to the sub-millimeter cores. 
Additionally, it is realized that observations at shorter wave- 
lengths can improve our understanding of the previously iden- 
tified star forming cores. This was the motivation to take a 
closer look in the mid-infrared (roughly 3-30 jum), where Young 
Stellar Objects (YSOs) are known to emit the bulk of their en- 
ergy. This paper describes the interpretation of Spitzer Space 
Telescope photometry of G48.65, focussing on those sources that 
show emission at 24yum, where the brightest, massive, embedded 
young stars are expected to be found. 

Section l2] describes the Spitzer data and their reduction. 
Section l3] deals with the classification of twenty sources based 
on the available data points in the near- and mid-infrared and 
sub-millimeter regimes. Conclusions and a discussion are pre- 
sented in Sect. |4] 
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Fig. 1. Twenty point sources visible in the 24 /^m MIPS image. 
Axes are in J2000 coordinates. The brightness scale is logarith- 
mic, ranging from 27 (black), through red and yellow, to 170 
MJy/sr (white). The "S" numbers correspond to the source IDs 
as listed in Table [T] The annuli around the brightest sources are 
artefacts of the instruments point spread function. 



2. Observations 

2. 1 . Spitzer photometry 

Mid-infrared images of IRDC G48.65 were recorded in October 
2004 by the Infrai-ed AiTay Camera (IRAC) and Multiband 
Imaging Photometer (MIPS) on board the Spitzer Space 
Telescope, as part of the proposal by Kraemer et al. (General 
Observing Proposal #3121). Photometric images were obtained 
in five broadband filters: all four bands of IRAC (centered 
around 3.6, 4.5, 5.8 and 8.0yL(m) and the first band of MIPS 

For the IRAC bands, the 12 second High Dynamic Range 
mode is used: five dithered images of 10.4 s and five of 0.4 s, 
adding up to a total integration time of 52-1-2 seconds per IRAC 
band per sky position. IRAC observes two adjacent 5' x 5' fields 
simultaneously, MIPS observes only one at a time. MIPS images 
are also dithered, a total of 44 images spread over the two fields 
add up to an integration time of 55 s per field. The resolution in 
band 1 of IRAC is limited by the pixel size of the CCD (~l'.'2); 
it is diffraction limited for the higher wavelengths, with the res- 
olution element ranging from r.'3 for IRAC 4.5 /zm to 6" for the 
24yumMIPS band (Rieke et al.,2004, .Fazio et al.|2004 1. 

The Spitzer Science Center (SSC) provides 'basic caUbrated 
data' in units of physical flux per unit solid angle. The basic 
calibrated data originate from pipeline version S 11 .4.0 in case 
of IRAC and version S 10.5.0 for the MIPS band. These im- 
ages are processed and combined into mosaics using scripts 
from the MOPEX software package ( Makovoz & Marleau 2005] 
Makovoz & Khan 2005). Data from the 0.4 s exposures is used 
to replace saturated pixels in the 10.4 s images. 



2.2. Point source extraction 

At Galactic coordinates (£, b) = (48?66, -0^30), IRDC G48.65 
lies close to the line of sight toward the giant molecular cloud 
(GMC) W51. Along a line of sight this close to the Galactic 
plane in general and the bright GMC in particular, there is signif- 
icant background radiation in the mid-infrared. The background 
varies with amplitudes comparable to some of the point source 
peaks and on size scales smaller than the point spread functions. 
This complicates (automated) point source extraction, especially 
at the longer wavelength bands at 8.0 //m and 24 jum. In the 
shorter wavelength bands, spatial variability of the background 
is less of an issue, but source confusion starts to play a role. 

The MIPS image (at 24yum) shows 20 distinct point source 
objects (see Fig.fTi. A subpixel position is determined for each 
of the 20 sources by fitting radial profiles in the 24yum image. 
Point spread function (PSF) photometry is performed at these 
positions, using the 'optimal photometry' routine ( |Naylor|1998[ l 
built into the STARLINK GAIA tool. The PSF itself is modeled 
on source S2, being bright and relatively isolated. The PSF is 
truncated at a radius of 6'.'1, the sky anulus is defined to be be- 
tween 12'.'3 and 15'.'3. An aperture correction factor of 1.7, ap- 
propriate for this radiu^ is applied. Sources S12 and S15 re- 
quire a smaller truncation radius (3'.'9) in order to prevent over- 
estimating the background sky level; in this case, an aperture 
correction factor of 2.5 is applied. The uncertainties listed for 
every lAfjxa. flux value in Table [Tlinclude sky pixel variance and 
the 4% uncertainty in the absolute calibration ( jEngelbracht et al. 
2007) . 



' |http : //ssc ■ spitzer . caltech . edu/mips/apercorr/ 
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The positions of the 24 jum sources are used as starting posi- 
tions for the photometry in the four IRAC bands. We apply the 
same method as described above, the only differences being: 

(i) the PSF is modeled on a bright isolated star in each image 
(S2 is not a suitable candidate at IRAC wavelengths); 

(ii) the PSF is truncated at 4'.'8 for IRAC 3.6 and 4.5 yum and at 
5" for IRAC 5.8 and 8.0 A^m; 

(iii) the relevant aperture correction is roughly 1.1; 

(iv) the absolute caUbration uncertainty is taken to be 5% ( |Reach| 
let al.|2005] l. 

The fluxes of the 20 MIPS sources, their counterparts in the 
IRAC bands, and the corresponding uncertainties are listed in 
Table [1] Sources S4, S17, S18 and S19 suff'er from bad pixel 
columns at the two shortest wavelengths; for these sources, the 
fluxes and uncertainties at 3.6 and 4.5 /im are taken from the 
GLIMPSE archive. At the positions where no counterpart was 
found for a MIPS source in a particular IRAC band, an upper 
limit to the flux value, indicated by a "<" symbol, is listed in 
Table [T] Since we want to avoid underestimating upper limits 
to the flux values in crowded, embedded and extincted regions, 
we conservatively set the upper limits around once or twice the 
lowest flux value that is detected in the same band. Here, the 
uncertainties are estimated at 50% of the flux limit value. 



2.3. Additional data: sub-millimeter, far- and near-infrared 

In addition to the newly obtained Spitzer Space Telescope mid- 
infrared photometry, we use sub-millimeter data points from 
lOrmel et al.K 2005 ) obtaine d with SCUBA at JCMT, MIPSGAL 
70 /zm far-infrared images (Carey et al. 2006|l, and archiva l near- 
infrared data from the 2MASS catalog ( [Skrutskie et al.||2006 i. 
See the corresponding columns in Table Ml for the objects that 
have counterparts in the near-infrared (2MASS) or in the sub- 
miUimeter (SCUBA). 

The SCUBA data point named PI is matched to S5 and S6 
and is discussed in Sect. |3.2.1| The positions of the peaks P2 
and EP as listed in Table 2 of |Ormel et al.| ( |2005 1> are incor- 
rect and have an offset of ~15" in RA with respect to the true 
peak positions in the 850 //m and 450 jiva images (C.W. Ormel, 
priv. comm., see also Fig.|2]l. The derived fluxes of the two peaks 
are nevertheless coiTect. The coiTect peak position for P2 is (RA 
= 19*'21™48'!0, Dec = +13°49'21") (±2", J2000); the correct 
peak position forEP is (19*'21'"44!9, -t-13°49'24") (±5", J2000). 
The mid-infrared source listed as S12 in Table[T]can therefore be 
matched to EP in the sub-millimeter, and S15 to P2. 

The third brightest 24fj.m object, S2, is also seen in the 
450 yum SCUBA image, albeit less bright than the other three 
peaks. Since it is not in the list of SCUBA peaks of |Ormel et al.| 
(2005 I, we estimate its 450 //m flux to be between 0.1 and 0.7 
times that of P2, i.e. 1.05 + 0.8 Jy (it is listed as 'P3' in Table[l). 

The MIPSGAL team recently released mosaics of the 70/im 
survey. We use the mosaics that overlap with our field of view to 
find one emission peak of roughly 7 Jy, coiTesponding to the S5 
and S6 mid-infrared sources (see Sect. |3.2.T] l. There are no other 
70 //m peaks, yielding a conservative upper limit of 5 Jy for the 
other sources in Table [T] 

Ten sources have counterparts in one or all of the near- 
infrared J, H and Ks 2MASS bands ( [Skrutskie eraT] [2006) 1, 
see Table [T] For the objects without 2MASS counterparts that 
show ambiguous model fitting results using only the Spitzer data 
points (S5, SIO, S16 and S20), 2MASS upper Hmits (magnitude 
15.8, 15.1 and 14.3 for J, H and K^ bands, respectively) are in- 
cluded as constraining parameters. 




Fig. 2. Composite image of IRDC G48.65. The color image 
is composed of images of the Spitzer passbands at 24 /urn 
(red), 8.0/im (green) and 4.5 yum (blue). Overlaid are the 
850jum SCUBA contours at levels of 0.08, 0.16, 0.24, 0.32 and 
0.40 Jy/beam. The axes annotate right ascension and declination 
(J2000). Notable MIPS point sources are labeled. 



3. Results 

3.1. Fitting data points to models of Young Stellar Objects 

The physical and chemical parameters and processes that play a 
role in specific YSO candidates are best explored by comparing 
the data points (flux values at various wavelengths) to models of 
YSOs. We use the models from Robita ille et al.[ ( [2006|l, who used 
a numerical Monte Carlo radiation transfer code d Whitney et al. 
|2003a|b[ [2004| l to calculate a set of 20000 models for YSOs. 
In this grid of models, 14 physical parameters characterize the 
three main components of a YSO: (i) the central star, (ii) a ro- 
tationally flattened infalling envelope with bipolar cavities, and 
(iii) a flared accretion disk. The effective number of parameters 
is 8 for models that have an envelope but no disk, and 10 for 
models that have a disk but no envelope. Apart from the 14 pa- 
rameters used as input for the radiative transfer code, 10 dif- 
ferent viewing angles are considered, ranging from edge-on to 
pole-on with respect to the axisymmetric geometry of all mod- 
eled systems. Every viewing angle results in a different Spectral 
Energy Distribution (SED) of the same model; the grid there- 
fore encompasses 10x20000 different SEDs. Various star for- 
mation 'modes' and stages, as well as environments, are taken 
into account by using a large spread in the input parameters. The 
code is developed bearing in mind the specific goal of applying 
the model grid to archives of Spitzer Space Telescope data, such 
as those of GLIMPSE (Benjamin et al.[[2003] l and MIPSGAL 
(Carey et al. 2006). This makes the model grid particularly use- 
ful for the interpretation of our data. 
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An SED fitting tool ( [Robitaille et al.||2007| ) is available on- 
lina^ Based on user-specified data points (flux and flux uncer- 
tainty in at least three passbands), the tool determines which of 
the models from the pre-computed grid ( [Robitaille et al.|2006| l 
fits best to the data. Upper limits in certain bands may also be 
specified, but do not contribute to A^, the number of data points. 
In addition, the user is asked to specify a fiducial distance range 
and a range in visual extinction. Ay. 

For every model in the grid, the specified measured fluxes 
(Fy{Ai)) are compared to convolved, scaled and extincted model 
fluxes (Py{Ai)), using x^ as a measure for the deviation of a cer- 
tain model from the specified data points: 



N 



1 '^ I 



<logio(F,(^,))) - logio(P,(A,)) 
(r«logio(F,(i,)))) 



(1) 



The o" factor in the denominator of Eq. ([TJ is the square root of 
the 'unbiased' flux variance in "^ log space, which is related to 
the relative flux uncertainty in each data point supplied by the 
user. The best fit is considered to be the SED that results in the 
minimum value of ;t'^ for the given set of data points: (;t^/^)best- 
An object can only be fitted to the model grid if A^ > 3, not 
counting upper limits. 

In addition to a comparison of a set of measured fluxes to the 
model YSO SEDs, the same data points are fitted to stellar pho- 
tosphere SEDs with effective temperatures ranging from 2000 to 
50 000 K. This feature is useful in assessing whether a partic- 
ular object might be explained by a normal stellar photosphere 
instead of a YSO model. 



3.2. Spectral Energy Distribution analysis 

A total of 17 of the 20 sources from Table [T] aided by sub- 
millimeter, MIPS 70 jum and 2MASS data, have A^ > 3 and 
can be analyzed using the models and the fitter from Robitaille 
[eraL] ( |2006,,2007j described in Sect.|TT] Three sources, SI, S3 
and SIO, do not have enough data points to use the fitting rou- 
tine. This analysis results in four suspected photospheres and 13 
YSOs, summarized in Table |2] 

In order to separate YSOs from (extincted) photospheres, the 
data points of the 17 sources are compared first to photosphere 
models. Visual inspection of all SEDs identifies the objects S4, 
S17, S18 and S19 (SEDs are shown in Fig.[3]l as photospheres. 
The x^ value for the best fitting photosphere is considerably 
lower than for the best fitting YSO model for all four sources. 
These sources are suspected to be foreground sources along the 
line of sight toward G48.65. None of the photosphere objects 
seem to be associated with the filament of the IRDC in the plane 
of the sky (see Fig.[T]l. Without any pre-conceived selection with 
respect to sky position, the other sources - the YSO candidates - 
have a location on the sky closer to the dark cloud, where we ex- 
pect star formation to occur. 

While AGB stars and even background galaxies are candi- 
dates to explain infrared excess in point sources, both types of 
objects are statistically unlikely to be present in the field of view 



considered in this paper. Using the model by Wainscoat et al. 
( |1992| l for the spatial distribution of AGB stars in the Galactic 
disk, a crude volume integration along the pyramid subtended 
by the lines of sight toward the 5' x 5' field yields an estimate 
of less than 1 AGB object in our field. A search through the 






Fig. 3. Best fitting photosphere model (solid lines) to the data 
points (filled circles) of S4, S17, S18 and S19. 

NAS A/IR\C extragalactic databasajresults in two known galax- 
ies within our field of view, neither of which coincide with any 
of the objects listed in Tablefl] Considering this, we assume here 
that all objects in our field that show infrared excess are in fact 
YSOs. 

The data points of the remaining 13 sources are compared to 
the YSO model grid of .RobitaiUe et al.| ( ,2006) (see Sect.[TT]i. We 
attempt to derive an evolutionary stage for each source, based on 
three key parameters of the fitted models: the rate of mass accre- 
tion from the circumstellar envelope (Mem) and the mass of the 
circumstellar disk (Mdisk), both relative to the mass of the cen- 
tral object. Mi,. We prefer the classification of YSOs into evo- 
lutionary "Stages" from Robitaille et al. ( 2006 1 over the older 



"Class O/I/II/IH" classification scheme, which is based solely on 



the slope of the mid-infrared SED (Adams et al.||1987 White 
et al^ |2007| l. The Stage classification provides a more physical 
basis for YSO classification: a source is assigned to a certain 
stage based on parameters of the model(s) that fit best to the 
data points. The threshold values for M^^ylM^, and M^isk/A^* 
are given in Table [3] We choose to make no distinction between 
"Stage 0" (with an SED resembling a 30 K graybody) and "Stage 
I" (a protostar with a large accreting envelope), since the en- 
ergy distributions of the two are difficult to distinguish in the 
mid-infrared; these two phases will be collectively referred to as 
"Stage I" in the rest of this paper. A Stage I object is dominated 
by the accreting envelope; the envelope accretion has diminished 
significantly in the consecutive Stage II, which shows signs of a 
massive circumstellar disk emerging from the envelope; in Stage 
III, the disk has become less and less massive and the direct ra- 
diation from the central star begins to dominate the observable 
flux. 

Tablel2]lists the 13 objects that are classified as YSOs. The 
spread in model parameters M*, Menv and Mdisk given in the ta- 
ble is obtained by considering all models that obey the criterion: 






+ 1, 



(2) 



with x^ IN as defined in Eq. (m. This criterion is chosen after 
visual inspection of the sets of models that fit the data points of 



http : //caravan . astro . wise . edu/protostars/ 



^ The NASA/IPAC Extragalactic Database (NED) is operated by the 
Jet Propulsion Laboratory, Cahfornia Institute of Technology, under 
contract with the National Aeronautics and Space Administration. 
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Fig. 4. YSO model SEDs that fit best to the data points (filled circles) of thirteen sources. This figure shows those sources listed in 
Table|2] Objects S2, S5, S6, S7, S9 and S12 are discussed in detail in Sect. [3XT] objects S8, Sll, S14 and S16 in Sect. [TZ2l and 
objects S13 and S20 in Sect. 3.2.3| In each panel, the black or colored curve is the best fitting SED (if it is colored, this indicates 
dependence on aperture size, see legend in bottom right frame). The lighter shaded curves (gray or pale colors) are all other model 
SEDs for which the criterion in Eq. (J2]i is obeyed. The dashed curve is the photosphere used as input for the radiative transfer code 
that calculated the best fitting YSO SED. 
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Table 2. Sources classified as YSO. 



(1) (2) 
ID assoc. 



(3) (4) (5) 

N minimal ;f^/A' M-,, 



(6) 

(Mo) 



1-3 X 10-* 
1-6 X 10-5 
5 X 10-'-2 X 
0-5 X 10^** 
2 X 10-'-3 X 
0-2 X lO-'* 
1-4 X lO-'* 
~ 10^' 
0-4 X 10^5 
2 X lO^" 
0-3 X 10-* 
~0 



(7) 

A^disk/M* 

(yr-') 



(8) 
Stage 



3 X 10-^-7 X 10"- 
1-5 X 10-- 
3 X 10-'^-5 X 10-- 
8 X 10-5-6 X lO-"* 
5 X 10-M X 10-^ 
5 X 10-«-7 X 10-- 
5 X 10-'*-2 X 10-' 

2 X 10-*-3 X 10-- 
-10-2 

3 X 10-'*-5 X 10-- 
3 X 10-''-5 X lO-"* 
3 X 10-M X 10-2 
-10-5-10-' 



S2 


yes 


7 


1.38 


S5 


yes 


8 


6.18 


S6 


yes 


7 


6.77 


S7 


yes 


4 


0.0875 


SB 


yes 


8 


1.99 


S9 


yes 


8 


0.126 


SU 


yes 


5 


0.61 


S12 


yes 


4 


0.58 


S13 


yes 


5 


11.9 


S14 


yes 


7 


0.030 


S15 


yes 


4 


0.515 


S16 


yes 


5 


1.99 


820 


no 


5 


18.6 



1-5 

1-2 

0.7-8 

0.2-6 

1.5-4 

0.3-6 

0.1-5 

1.5-6 

-2 

0.5-4 

3-4 

2-5 

-4-6 



10- 



10- 



I 

I 

I 

I 

II 

I or (II) 

I or II or (III) 

I 

-I 

(I) or II or (III) 

I 

I or II 

-II 



Object parameters derived from fitting the data points to YSO SEDs (see Sec. |3.2| (. Column (1) lists the ID of the object, corresponding to 
the entry in Table [T| Column (2) indicates whether the source is thought to be spatially associated with the IRDC filament. The number of 
available data points, A', is listed in column (3). The x^ per data point, x^/^' is given in column (4). Columns (5), (6) and (7) list the fitted 
stellar mass, envelope accretion rate and disk mass, respectively. Note that the latter two are given in terms of M-^, not Mq. Finally, column 
(8) lists the Stage in which the sources fall according to the values in columns (6) and (7) and the Stage definition in Table [3] If multiple 
possibilities are given for the Stage of a source, the less likely one(s) is (are) presented in brackets. A '-' prefix in column (5), (6), (7) or (8) 
indicates a high uncertainty in the derived parameter, resulting either from a high value of the minimal ^^/yy and/or a large spread in fitted 
parameter values. 



Table 3. YSO evolutionary stages as defined by Robitaille et al. 
( |2006l l. 



Stage components 



Md,sk/M* 



Me„v/M* 



I accreting envelope any > 10- yr- 

II disk, remains of envelope > 10-* < 10-'' yr- 

III optically thin disk < lO-'' < 10-'' yr- 



each of the objects: it should not be so stringent that only one 
model fits a source, and on the other hand, it should not be so in- 
discriminate that it allows the majority of the models in the grid 
to be a "good fit". It is subsequently kept fixed for all sources. If 
Table |2] lists only one value for a certain parameter, the value is 
implied to be known to an accuracy indicated by the number of 
significant digits. We refrain from defining a proper confidence 
level for each derived parameter value, since the 14 dimensional 
space of physical parameters is too sparsely sampled to derive 
any formal confidence intervals ( ,Robitaille et al.,2007) . The ob- 
jects classified as YSOs are discussed in detail below, starting 
with the Stage I objects. 

3.2.1. Stage I: S2, S5, S6, S7, S9, S12 and S15 



The sub-millimeter core "PI", identified by |Ormel et al.| ( [2005| l, 
is separated into two distinct sources by MIPS and IRAC (see 
Fig.l2]l. What Ormel et al. (2005 1 modeled as one emission core, 
now turns out to be composed of at least two individual YSOs: 
S6 and S5 in this paper. The 24 fim peaks of these sources are 
separated by roughly 12", about 0.15 pc tangential distance at 
2.5 kpc distance. It is crudely assumed that -50-80% of the 450 
and -50-90% of the 850 //m SCUBA flux of PI originates from 
S6, the brightest of the two cores in 24 fim; the remainder of 
the PI sub-millimeter flux is attributed to S5, with the same un- 
certainty margins. The large uncertainties are adopted mainly to 
reflect the uncertainty in the assumption of the relative contribu- 
tions of S5 and S6 to the total PI flux. In addition, the MIPSGAL 



70 //m mosaic image j Carey et al.||2006 l that covers G48.65 is 
found to contain an emission peak at the position of PI (S5 and 
S6). Its flux is measured using the same method as for the IRAC 
and MIPS 24 jum bands (see Sect. 2.2 1. Its flux is found to be 
6.9±0.3 Jy at 70 yum. Like for the sub-millimeter flux described 
above, it is uncertain how much of this flux is contributed by 
each individual mid-infrared source (S5 and S6 fall inside one 
MIPS 70 fixn beam). The 70 fixn flux is therefore attributed to S5 
and S6 in the same proportion as the 450 |/m flux, i.e. -50-80% 
to S6 and -20-50% to S5. 

Both S6 and S5 (see Fig.|4]for the SED fits) are fit by models 
with high envelope accretion rates. Figure [5] shows the spread 
in Menv-A^disk spacc of models fitted to S6; the diagram looks 
similar for S5. Both sources are therefore classified as Stage I 
objects. The total luminosity according to the best-fitting models 
is -10^-10^ Lq for S6 and -10^ L© for S5. The sum is consistent 
with t he 1-cr limit of the total luminosity that was modeled for 
PI by 



Ormel et al. 



(2005):330|g^o"'L, 



The sources S15 (sub-millimeter peak EP in |Ormel et al.| 



2005 1 and S 12 (P2) are only detected in MIPS 24 fim, in the sub- 



millimeter, very faintly in the IRAC 4.5;um band, and not in the 
other IRAC bands. Like S5 and S6, they are well fit by models 
with high values of Menv, exceeding 10-^ Mo/yr. FigureHshows 
the SED fits for S15 and S12. The inferred stellar masses for 
both objects are between 1 and 6 Mq, the total luminosities are 
-10^ Lq each. Although S12 (EP) appeared slightly less lumi- 
nous from sub-millimeter continuum modeling, the rates of en- 
ergy production inferred here for S 12 and S 15 are largely consis- 
tent with the modeled l uminosity values for EP and P2 in Table 



3 of Ormel et al. 



(|2005|, 300!'™" and 19_i7 



110 



Lq, respectively. 



Model fitting for the S2 data points, including the 450 //m 
data point (see Sect. 2.3 i, classifies it as a 1-5 Mq Stage I object, 
see Table|2] S7 and S9, on the northern side of the cloud filament, 
are also classified as a Stage I objects in Table |2] A minority of 
the well fitting models lie in the Stage II region for both these 
objects. Nevertheless, it is deemed more likely that they are in 
Stage I. 
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Fig. 5. Ranges of Menv (vertical axis) and Mdisk (horizontal axis) 
of models that fit the data of S6. The gray levels show the distri- 
bution of the models in the grid, where darker color indicates a 
higher density of models. The points show the loci of all models 
which obey the criterion in Eq. (|2]), i.e. the model SEDs which 
are plotted in Fig. HI 



Summarizing, all mid-infrared objects that show counter- 
parts at 850 yum and/or 450 //m represent the very earliest stage 
of star formation. However, S7 and S9 - with no data points 
longward of 30 yum - show that a sub-millimeter counterpart is 
not a required condition for a source to be classified as Stage I. 



3.2.2. Stagell:S8, S11,S14andS16 

S8 is found to match models with a stellar mass of 1 .5^ Mg, and 
no accretion from the envelope. The absence of envelope accre- 
tion, shown by the distribution in Menv-A^disk space in Fig. 6] 
rules out the option that SB is an embedded Stage I source. The 
relative disk mass is high enough to rale out Stage III and S8 is 
therefore classified as a Stage II object. 

Objects Sll, S14 and S16 are all intrinsically faint in the 
IRAC bands. It is evident from the SED plots of these objects in 
Fig. fflthat model fitting therefore benefits less from upper limits 
from MIPS 70 fim and SCUBA 450 yum and 850 i^m. The lack 
of data points and stringent upper limits at far-infrared and sub- 
millimeter wavelengths results in an ambiguity: the presence of 
the bump longward of ~30 yum, a characteristic feature in SEDs 
of Stage I sources, cannot be ruled out. This is reflected by the 
envelope accretion rate in column (6) of Table l2] which ranges 
from to well above the Menv - lO^^M^^yr threshold. For 
S14, the majority of models that obey the criterion in Eq. l2]use 
an envelope accretion rate below this threshold, and it is there- 
fore classified to be most likely a Stage II object, with a mass 
between 0.5 and 4 Mq. Sll and SI 6, however, show a partic- 
ularly large spread of models covering both Stage I and Stage 
II parameters, which is why their evolutionary stage in Table |2] 
is listed as "Stage I or 11". The models fitted to Sll and S16 
have stellar masses ranging from 0.1 to 5 and from 2 to 5 Mq, 
respectively. 
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3.2.3. S13andS20 

The criterion in Eq. l2]leaves only one YSO model for S13, and 
even this is a relatively poor fit (see Fig. |4] and column (4) of 
Table l2]). Hence, no firm conclusion can be drawn regarding any 
of the key physical parameters such as stellar mass, disk mass, 
envelope accretion rate, and therefore its evolutionary stage. 

Finally, S20 is the only source classified as YSO which is 
not spatially (in the plane of the sky) associated with the IRDC. 
Its key parameters are poorly constrained due to the high devi- 
ation between the best-fitting model and the data points of S20 
(see Fig.|4|. Therefore, the derived parameters and evolutionary 
classification of this objects should be considered with care. 



3.3. Color Analysis of S1 and S3 

Since SI, S3 and SIO have less than 3 data points, it is noted in 
Sect. l3.1l that model SEDs cannot be fit to these sources. Source 
SIO is overwhelmed by diffuse background radiation. We are 
therefore reluctant to draw conclusions about the physical prop- 
erties of this object based on its two data points. 

Sources SI and S3 on the other hand, although they only 
have a flux measurement at 24 fj.m, can be examined purely by 
their [8.0] - [24] color. Using the MIPS zero point magnitudes 
from the SSC website^ SI and S3 have a [24] magnitude of 6.1 
and 5.8, respectively. Using the 8.0 yum upper limits from Table [T] 
and the IRAC magnitude system defined in Reach et al. ( 2005|l, 
the [8.0] magnitude is at least 11.3 for each source. This puts a 
lower limit to their [8.0] - [24] colors: > 5.2 for SI and > 5.5 for 



S3. Figure 23 from Robitaille et al. ( 2006 1 shows that practically 
only Stage I YSO models exhibit an [8.0] - [24] color this red. In 
conclusion, if an evolutionary phase is to be attributed to sources 
S 1 and S3 based on the above argument, it would be Stage I. 



|http : //ssc ■ spitzer . caltech . edu/mips/calib/ 
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4. Conclusions and Discussion 

4.1. Conclusions 

Of the 20 sources near IRDC G48.65 that are visible at 24 //m, a 
total of 13 are classified as YSOs, seven of which are found to 
be in Stage I, two are in Stage II, two more are in either Stage 
I or Stage II, and two are uncertain. While all reliably classified 
YSOs lie along the IRDC extinction filament, the four sources 
matched to photosphere models all lie away from the filament. 
Each YSO generally resides within a projected distance of <1 pc 
from another YSO. We conclude that G48.65 is an example of a 
dark cloud environment forming a group of stars. 

The most important modeled properties of the 24 /im sources 
that are classified as YSOs are summarized in Table |2] The stel- 
lar masses of the modeled objects range from slightly sub-solar 
to ~8 Mq. The inferred evolutionary Stages are predominantly 
early phases, ranging from early Stage I to late Stage II. This is 
believed to be caused partly by the selection effect discussed in 
Sect. 4.4 Therefore, the presented range of evolutionary phases 



should not be interpreted as a full picture of all star forming ac- 
tivity in this IRDC. 

The emission peak PI identified by Ormel et al.| ( 2005 i is 
resolved by Spitzer into two distinct emission cores: S6 and S5. 
Both objects are found to be in the earliest Stage of star forma- 
tion. Stellar masses are not well constrained for these sources 
but are unlikely to be in excess of 8 Mq. The summed total lumi- 
nosity of the best-fit models for S5 and S6 is ~10^-10^ Lq, con- 
sistent with results from modeling by |Ormel et al.| ( |2005| based 
purely on sub-mm observations. The total luminosity of the best 
fitting model for objects S12 (EP) and S15 (P2) are consistent 
with values found by [Ormel et al. (2005[l. 



4.2. Robustness of Results 

It is important to assess a YSO based on data points in a wave- 
length regime as extended as possible. While a detection of a 
certain YSO at sub-millimeter wavelengths seems to be a strong 
indication that it is in Stage I (S5, S6, S12, S15), a detection 
by 2MASS in the near-infrared does not rule out a Stage I clas- 
sification, as sources S2 and S7 show. It may therefore lead to 
erroneous classification of YSOs if one bases it either on only 
near-infrared or on only far-infrared/sub-millimeter data. 

The availability of data points in the 3-30 /vm regime is in 
some cases enough to constrain the evolutionary stage of a YSO. 
Source S8 (see Fig. HI for the SED fit) is an example in which 
the degeneracy of the model fitting is decreased considerably 
by adding 2MASS data points (see also Sect. 3.3 and Fig. 3 of 
[Robitaille et al., ( ,2007j ). 

4.3. Continuous Low-mass Star Formation 

The Stage I sources in particular show large uncertainties in stel- 
lar mass (see Table l2]). This is not surprising considering that, 
for objects that emit mostly through re-radiation from the ac- 
creting envelope, it is of course inherently diflicult to model pa- 
rameters of the obscured central heating source. However, the 
OH maser survey by |Pandian et a l. (2007 ) indicates the absence 
of OH maser emission in the direction of IRDC G48.65. Since 
OH masers are commonly associated with high-mass star forma- 
tion (>8 Mq), this agrees with the fact that we find no massive 
YSOs in G48.65. Considering the above, we are confident that 
any YSO associated to G48.65 is unlikely to harbor a central star 
more massive than 8 Mq. 



Based on the range of stellar masses and evolutionary Stages 
of the YSOs, we conclude that G48.65 has been an active low- 
mass star forming region for at least ~10* years (late Stage II) up 
to as recent as ~10"* years ago (early Stage I) and is likely to be 
in the process of forming still younger stars. Assuming that the 
IRDC is indeed the birth place of these young stars, the cloud 
itself must have been stable over a time scale of at least ~10* 
years. This is consistent with modern views on lifetimes of <10* 
years for molecular clouds (as reviewed in Mac Low & Klessen 
2004') and with lifetimes of <10^ years for dense cluster forming 
clouds (see review by Larson|2003[l. 



4.4. Mid-infrared Selection Effect 

The starting point of the study in this paper - the set of point 
sources in the MIPS 24jum image - introduces an observational 
bias toward the younger evolutionary stages. Early stage YSOs 
are simply brighter in the mid-infrared, since the total mid- 
infrared luminosity of a YSO generally decreases with time. 
Especially in a fairly distant region such as IRDC G48.65, at 
~2.5 kpc, one must keep in mind that later stages and intrinsi- 
cally fainter sources are less likely to be detected. This study in 
particular is biased toward the earlier stages, since it focusses on 
the 24yum objects. 

It is primarily the relatively cold envelope component, dom- 
inant in Stage I objects, that gives a YSO SED its large "bump" 
longward of 20/vm. Hence, an analysis of objects near the IRDC 
that are only visible at wavelengths shortward of 8 fj.ni is ex- 
pected to lead to the identification of more evolved (Stages II 
and III) and less massive YSOs. The IRAC 8 fim image, for 
example, shows at least a factor 3 more point sources than the 
24 fj.m image in the same field of view. Such a study will be es- 
sential in getting a complete census of the star formation activity 
in G48.65, both in terms of time scales and in terms of a stellar 
mass function. 

In principle, the set of young stars of various masses pre- 
sented in this paper could be used to construct an initial mass 
function (IMF). Extrapolation of such an IMF could provide in- 
sight into the absence of massive stars (>8 Mq). However, we 
refrain from presenting such an exercise here. The first reason 
for this is the incompleteness of our sample discussed above; a 
second reason is that the inferred stellar mass of each individual 
YSO (see Col. (5) of Tablel2]i is rather poorly constrained. 

4.5. Outlook 

In the near future, the 2MASS point source catalog (l-3;um) 
and the GLIMPSE catalog (3-9.5 /im) will be complemented 
by the MIPSGAL point source list of the Galactic plane at 24 
and 70 //m. Moreover, the UKIRT Infrared Deep Sky Survey 
( [Lawrence et al.|2007| l is already partly available to the general 
community and will eventually cover half of the Galactic plane 
in the /, H and K bands, superseding the sensitivity of 2MASS 
by roughly 3 magnitudes. This will open up the possibility to in- 
vestigate YSOs in hundreds of other IRDCs in the Galaxy using 
an approach similar to the one taken in this paper. This approach 
has shown that source-by-source investigation of emission cores 
with l-30yum data can lead to the identification of young stars 
in various phases of evolution. Perhaps positions that show sub- 
millimeter emission are the most interesting cases to start with 
if one is attempting to identify very early stages of star forma- 
tion. However, even with the combination of data from 2MASS, 
GLIMPSE and MIPSGAL, covering just the 1-70 //m regime, it 
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is possible to systematically study star formation activity across 
the Galactic plane. IRDCs may eventually be divided into groups 
according to whether they show star formation activity, and if 
they do, whether they form low-mass stars or high-mass stars. 

The Stage 1 sources in Table l2] are expected to show associ- 
ated outflows and shocks, traced by e.g. SiO. Detection of such 
tracers could confirm the star forming activity and reveal kine- 
matic and spatial structure of YSOs. In addition, measurements 
by the VISIR instrument (8-13jum) at VLT will be able to re- 
solve more detailed structure (<10^ AU at 2.5 kpc) of the sources 
and show envelope and disk components. It may even show that 
objects that appear as one source in the Spitzer images are in fact 
multiple YSOs. Source S5, for example, appears slightly elon- 
gated in the Spitzer images. In fact, some of the poorer model 
fits (see Col. (4) in Table l2|i may be explained by the fact that 
the models do not account for multiple young star-disk-envelope 
systems inside one telescope beam. 

The high spatial and spectral resolution and new wavelength 
regimes of future observatories, such as ALMA, the Herschel 



Space Observatory (specifically the Hi-GAL survey, Noriega- 
[Crespo & MoUnari|2008) l and the James Webb Space Telescope, 
wiir enable more detailed studies of individual objects associated 
with G48.65 and other (distant) IRDCs. 

Acknowledgements. The authors are grateful to Thomas Robitaille for his help 
with the use of the YSO models, and for making adjustments to the fitting tool 
upon request. Furthermore, the authors wish to thank Floris van der Tak and 
Marco Spaans for their support and suggestions that were essential for the in- 
terpretation of the results and the completion of the manuscript. We thank Chiis 
Ormel and Wilfred Frieswijk for additional constructive discussions, comments 
and suggestions that helped to improve this paper We acknowledge the efforts 
of Kathleen Kraemer for preparing the observing proposal for the Spitzer data 
on which this paper is based and for her comments on the manuscript. We thank 
the anonymous referee whose comments and suggestions helped to improve the 
content as well as the clarity and presentation of this paper 

This work is based on observations obtained with the Spitzer Space 
Telescope, which is operated by the Jet Propulsion Laboratory, California 
Institute of Technology, under NASA contract 1407. In addition, this publica- 
tion makes use of data products from the Two Micron All Sky Survey, which 
is a joint project of the University of Massachusetts and the Infrared Processing 
and Analysis Center/California Institute of Technology, funded by the National 
Aeronautics and Space Administration and the National Science Foundation. 



Ormel, C. W, Shipman, R. F, Ossenkopf, V., & Helmich, R R 2005. A&A, 439, 

613 
Pandian, J. D., Goldsmith, P F, & Deshpande, A. A. 2007, ApJ, 656, 255 
Perault, M., Omont, A., Simon, G., et al. 1996, A&A, 315, L165 
Rathbome, J. M., Jackson, J. M., Chambers, E. T, et al. 2005, ApJ, 630, L181 
Rathbome, J. M., Jackson, J. M.. & Simon, R. 2006, ApJ, 641, 389 
Reach, W. T, Megeath, S. T, Cohen, M., et al. 2005, PASP, 1 17, 978 
Rieke, G. H., Young, E. T, Engelbracht, C. W., et al. 2004, ApJS, 154, 25 
Robitaille, T. P, Whitney, B. A.. Indebetouw, R., & Wood, K. 2007, ApJS, 169, 

328 
Robitaille, T. P., Whitney, B. A., Indebetouw, R., Wood, K., & Denzmore, P. 

2006, ApJS, 167,256 
Shipman, R. F, Frieswijk, W., & Helmich, F P 2003, in ASP Conf. Ser 287: 

Galactic Star Formation Across the Stellai' Mass Spectrum, ed. J. M. De 

Buizer & N. S. van der Bliek, 252-256 
Simon, R., Jackson, J. M., Rathbome, J. M., & Chambers, E. T 2006a, ApJ, 639, 

227 
Simon, R., Rathbome, J. M., Shah, R. Y, Jackson, J. M., & Chambers, E. T 

2006b, ApJ, 653, 1325 
Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1 163 
Wainscoat, R. J., Cohen, M., Volk, K., Walker, H. J., & Schwartz, D. E. 1992, 

ApJS, 83, 111 
White, R. J., Greene, T. P, Doppmann, G. W, Covey, K. R., & Hillenbrand, 

L. A. 2007, in Protostars and Planets V, ed. B. Reipurth, D. Jewitt, & K. Keil, 

117-132 
Whitney, B. A., Indebetouw, R., Bjorkman, J. E., & Wood. K. 2004, ApJ, 617, 

1177 
Whitney, B. A., Wood, K., Bjorkman. J. E., & Cohen, M. 2003a, ApJ, 598, 1079 
Whitney, B. A., Wood, K., Bjorkman, J. E., & WolS, M. J. 2003b, ApJ, 591, 

1049 



References 

Adams, F C, Lada, C. J., & Shu, F H. 1987, ApJ, 312, 788 

Benjamin, R. A., Churchwell, E., Babler, B. L., et al. 2003, PASP 115, 953 

Beuther, H. & Steinacker, J. 2007, ApJ, 656, L85 

Carey, S. J., Egan, M. P., Kuchar, T. A., et al. 2000, in Bulletin of the American 

Astronomical Society, 516 — h 
Carey, S. J., Noriega-Crespo, A., Mizuno, D. R., et al. 2006, in Bulletin 

of the American Astronomical Society, Vol. 38, Bulletin of the American 

Astronomical Society, 1023 — t- 
Egan, M. P, Shipman, R. F, Price, S. D., et al. 1998, ApJ, 494, L199+ 
Engelbracht, C. W., Blaylock, M., Su, K. Y. L., et al. 2007, PASR 1 19, 994 
Fazio, G. G., Hora, J. L., Allen, L. E., et al. 2004, ApJS, 154, 10 
Frieswijk, W. W. F, Spaans, M., Shipman, R. F, Teyssier, D., & Hily-Blant, P. 

2007, A&A, 475, 263 
Indebetouw, R., Robitaille, T R, Whitney, B. A., et al. 2007, ApJ, 666, 321 
Larson, R. B. 2003, Reports of Progress in Physics, 66, 1651 
Lawrence, A., Warren, S. J., Almaini, O., et al. 2007, MNRAS, 379, 1599 
Mac Low, M.-M. & Klessen, R. S. 2004, Reviews of Modem Physics, 76, 125 
Makovoz, D. & Khan, I. 2005, in ASP Conf. Ser 347: Astronomical Data 

Analysis Software and Systems XIV, ed. P. Shopbell, M. Britton, & R. Ebert, 

81-+ 
Makovoz, D. & Marleau, R R. 2005, PASP 117,1113 
Menten, K. M., Pillai, T, & Wyrowski, F 2005, in lAU Symposium, Vol. 227, 

Massive Star Birth: A Crossroads of Astrophysics, ed. R. Cesaroni, M. Felli, 

E. Churchwell, & M. Walmsley, 23-34 
Naylor, T. 1998, MNRAS, 296, 339 
Noriega-Crespo. A. & Molinari, S. 2008, in American Astronomical Society 

Meeting Abstracts, Vol. 212, American Astronomical Society Meeting 

Abstracts, 96.04-+ 



